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Global-scale  interactions  between  the  tropics  and  extratropics  are  investigated  using  a  version  of  the  I  S 
Navy's  global  operational  numerical  weather  prediction  model  The  primary  goals  of  this  study  are  I  )  to  dem¬ 
onstrate  the  importance  of  atmospheric  tcleconnections  for  medium-range  numerical  weather  prediction  and 
2 1  to  analyze  the  evolution  and  dynamic  structure  of  the  response  in  a  sophisticated  numerical  forecast  model 
The  model  normal  modes  are  used  as  the  principal  diagnostic  tool  for  analyzing  the  response  to  sea  surface- 
temperature  anomalies  in  the  tropical  Pacific.  By  monitoring  the  energy  growth  in  the  dominant  horizontal 
and  vertical  modes  and  comparing  these  with  conventional  difference-field  diagnostics,  it  is  shown  that  the 
character  of  the  long-term  response  is  well  established  within  one  to  two  weeks  after  the  heating  anomaly  is 
introduced.  The  growth  rates  and  structures  of  these  modes  provide  insights  into  the  dynamic  processes  that 
control  the  model  response.  In  the  tropics,  enhanced  convection  is  clearly  the  dominant  forcing  mechanism  for 
these  modes.  In  the  extratropics,  a  more  complicat'd  picture  arises  in  which  both  meridionally  propagating 
energy  and  in  situ  instabilities  in  the  ambient  flow  appear  to  be  important  mechanisms  for  producing  the 
observed  wave  patterns.  The  results  clearly  demonstrate  that  tropical  forcing  can  have  a  significant  global  impact 
on  time  scales  relevant  to  medium-range  numerical  weather  prediction. 

In  Part  II  of  this  study,  the  normal-mode  diagnostic  approach  is  extended  by  enveloping  a  technique  for 
partitioning  the  modes  according  to  their  latitudinal  variances  in  order  to  examine  the  tropical  and  cxlratropical 
responses  in  further  detail.  It  is  shown  that  the  modes  are  a  powerful  and  flexible  tool  for  diagnosing  the  behavior 
of  a  complicated  model. 


1.  Introduction 

There  is  no  doubt  that  tropical  forcing  can  have  a 
significant  impact  on  the  global  circulation  at  midlat¬ 
itudes  (e.g.,  Bjerknes  1969:  Horel  and  Wallace  1981: 
Wallace  and  Gutzler  1981;  Webster  1981;  Hoskins  and 
Karoly  1 98 1 ;  Geisler  et  al.  1985;  Pitcher  et  al.  1988). 
Such  interactions  between  the  tropics  and  extratropics, 
commonly  referred  to  as  teleconnections,  are  often  as¬ 
sociated  with  anomalous  variations  in  tropical  sea  sur¬ 
face  temperature  (SST),  such  as  those  that  occur  in 
conjunction  with  the  El  Nino-Southern  Oscillation 
(ENSO)  phenomenon  (Julian  and  Chervin  1978). 
There  is  now  growing  evidence,  however,  that,  in  ad¬ 
dition  to  changes  in  the  global-scale  standing  wave  pat¬ 
tern,  tropical  forcing  may  have  a  crucial  impact  on  the 
evolution  of  atmospheric  phenomena  on  a  wide  range 
of  spatial  and  temporal  scales.  For  example,  studies  by 
Kok  and  Opsteegh  ( 1985)  and  Mechoso  et  al.  ( 1987 ) 
indicate  that  stnaller-scale  transient  eddies  may  also 
be  affected  by  tropical  forcing  anomalies.  In  addition, 
it  has  been  shown  by  Tiedtke  ( 1984)  and  by  Donner 
et  al.  ( 1982)  that  the  response  of  general  circulation 
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models  (GCMs)  may  be  quite  sensitive  to  even  minor 
variations  in  model-predicted  diahatic  forcing  rates. 
These  and  other  results  suggest  that  interactions  be¬ 
tween  the  tropics  and  midlatitudes  occur  regularly  and 
are  an  important  component  of  the  dynamics  of  the 
global  circulation.  Where  numerical  weather  prediction 
(NWP)  is  concerned,  these  results  suggest  that  certain 
systematic  forecast  errors  at  midlatitudes  may  in  fact 
be  linked  directly  to  systematic  errors  in  model-pre¬ 
dicted  tropical  forcing. 

At  the  present  time,  much  of  what  is  known  about 
the  structure  and  dynamics  of  Vropical-midlatitude  in¬ 
teractions  has  been  learned  from  studies  involving  rel¬ 
atively  simple  models  whose  behavior  can  be  easily 
interpreted.  These  studies  have  shown  that  perturba¬ 
tions  in  the  tropical  atmosphere  caused  by  diahatic 
forcing  are  separable  into  two  basic  types  of  motions 
according  to  their  vertical  structures  (Lim  and  Chang 
1 983).  The  first  type  tends  to  be  trapped  at  low  latitudes 
and  is  characterized  by  cellular  motion  in  which  the 
large-scale  flow  is  of  the  opposite  sign  at  upper  and 
lower  levels.  The  second  type  is  characterized  by  large- 
scale  flow  that  has  the  same  sign  at  all  levels  and  can 
propagate  meridionally  from  the  tropics  to  high  lati¬ 
tudes.  The  latter  may  produce  the  global  teleconnection 
patterns  [ e.g.,  the  Pacific-North  American  ( PNA )  pat- 


4 


2898' 


MONTHLY  Wf  A  I  II!  R  R  1  V  I  I  W 


V.  H  l  Ml  l  ?<t 


tern,  Horel  and  Wallace  1981;  Wallace  and  Gutzler 
1981  ]  often  associated  with  the  atmospheric  response 
to  tropical  SST  anomalies. 

While  idealized  models  have  successfully  reproduced 
certain  qualitative  aspects  of  the  observed  response, 
much  remains  unknown  about  the  dynamics  of  this 
response  in  the  atmosphere,  which  may  be  complicated 
by  processes  such  as  waveguiding  by  subtropical  jets 
and  instabilities  in  the  ambient  flow.  Although  there 
now  exists  a  considerable  amount  of  data  documenting 
the  response  of  sophisticated  GCMs  to  anomalous 
tropical  forcing  (e.g.,  Keshavamurty  1982;  Blackmon 
et  al.  1983;  Shukla  and  Wallace  1983;  Palmer  and 
Mansfield  1986;  Mechoso  et  al.  1987;  Pitcher  et  al. 
1988).  little  has  been  presented  about  the  details  of 
the  modal  structure  and  the  dynamics  of  the  response 
in  these  models  or  about  the  extent  to  which  these 
characteristics  agree  with  results  from  idealized  mod¬ 
eling  studies.  In  addition,  most  of  these  studies  have 
examined  the  long-term  ( seasonal  to  interannual )  im¬ 
pact  of  tropical  forcing,  focusing  little  attention  on  time 
scales  relevant  to  NWP. 

In  this  study,  the  normal  modes  of  a  sophisticated 
NW  P  model  are  utilized  to  investigate  tropical-mid¬ 
latitude  interactions.  The  experimental  approach  is  to 
analyze  the  model  response  to  imposed  SST  anomalies 
in  the  tropical  Pacific.  The  main  objectives  are 

1 )  to  examine  the  medium-range  global  impact  of 
tropical  forcing; 

2 )  to  investigate  the  evolution  and  dynamic  struc¬ 
ture  of  this  response  in  a  realistic  system. 

A  secondary  objective  is  to  demonstrate  how  the  nor¬ 
mal  modes  may  be  used  effectively  to  diagnose  the  be¬ 
havior  of  complicated  models.  These  objectives  differ 
from  those  of  many  previous  GCM  studies  of  atmo¬ 
spheric  teleconnections,  which  have  provided  more  of 
a  climatological  overview  of  how  the  atmosphere  re¬ 
sponds  to  tropical  heating  anomalies. 

The  model  used  in  this  study  is  a  version  of  the  fore¬ 
cast  component  of  the  Navy  Operational  Global  At¬ 
mospheric  Prediction  System  (NOGAPS),  described 
by  Hogan  and  Rosmond  ( 1991 ).  The  results  are  re¬ 
ported  in  two  parts,  this  paper  being  the  first  part.  The 
analysis  in  this  paper  is  conducted  from  a  global  per¬ 
spective,  allowing  exploitation  of  the  simple  form  of 
the  normal-mode  energy  equations.  In  the  following 
section,  the  model  and  its  normal  modes  are  briefly 
described.  In  section  3,  we  describe  the  numerical  ex¬ 
periments  that  provide  the  data  for  this  study.  In  section 

4,  the  time-averaged  model  response  is  examined  using 
conventional  difference-field  diagnostics.  This  provides 
a  basis  for  comparison  with  previous  studies.  In  section 

5,  we  begin  the  normal-mode  analysis  of  the  global 
response.  The  evolution  of  ihe  response  is  analyzed  in 
terms  of  the  energy  of  the  modes.  In  section  6,  the 
dominant  horizontal  structures  of  the  response  are 
identified  and  their  evolution  is  examined.  Finally,  in 


section  7.  some  conclusions  and  speculations  based  on 
this  work  are  presented. 

In  Part  II  ( Gelaro  1 992;  hereafter  referred  to  as  Gil ). 
we  extend  the  analysis  in  this  paper  by  investigating 
the  tropical  and  extratropical  components  of  the  re¬ 
sponse  in  detail.  A  new  technique  is  developed  whereby 
the  normal  modes  are  partitioned  according  to  their 
latitudinal  variances,  allowing  us  to  analyze  the  evo¬ 
lution  and  structure  of  the  response  in  each  region  and 
to  make  inferences  about  the  dynamic  mechanisms 
that  produce  them. 


2.  Description  of  the  model  and  normal  modes 

The  model  used  in  this  study  is  version  3.0  of  NO¬ 
GAPS,  which  served  as  the  U.S.  Navy's  operational 
global  NWP  model  betwee  i  January  1988  and  March 
1989.  NOGAPS  is  a  complete  real-time  forecast  system 
consisting  of  an  optimal-interpolation  (Ol)  analysis 
scheme  ( Barker  et  al.  1989 ),  a  nonlinear  normal-mode 
initialization  (NNMI)  scheme,  and  a  global  forecast 
model  ( Hogan  et  al.  1991 ).  The  forecast  model  is  spec¬ 
tral  in  both  horizontal  dimensions,  with  a  triangular 
truncation  at  wavenumber  47  (T47  truncation ).  which 
corresponds  to  a  transform  grid  with  1 44  equally  spaced 
longitudinal  points  and  72  Gaussian  latitudes  ( roughly. 
2.5°  horizontal  resolution).  It  has  18  vertical  levels 
defined  in  terms  of  er  coordinates.  The  time-integration 
scheme  is  a  centered  semi-implicit  scheme  (Kwizak 
and  Robert  1971 )  with  a  Robert  ( 1966)  time  filter.  A 
biharmonic  (V4)  spatial  filter  is  employed  to  smooth 
extremely  small-scale  features. 

NOGAPS  contains  an  extensive  array  of  parame- 


Table:  I.  Values  of  7)  calculated  from  a  model  control  run.  and 
the  corresponding  set  of  equivalent  depths  h,.  Note  that  /  is  used  as 
a  model-level  index  for  a  and  f.  but  refers  to  the  ordering  of  the 
vertical  modes  in  the  case  of  h. 


1 

a 

f,<K) 

Mm) 

1 

0.008 

229.25 

9669.53 

2 

0.028 

209.14 

3224.16 

3 

0.049 

198.15 

801.71 

4 

0.092 

199.12 

308.66 

5 

0.136 

214.05 

130.55 

6 

0.187 

217.82 

63.85 

7 

0.247 

227.33 

40.60 

8 

0.315 

236.96 

21  30 

9 

0.393 

247.70 

11.48 

to 

0.482 

258.23 

663 

11 

0.575 

266.92 

3.99 

12 

0.669 

273.94 

2.36 

13 

0.759 

279.79 

1.45 

14 

0.839 

283.94 

0.75 

15 

0.904 

287.66 

0.40 

16 

0.950 

290.40 

0.18 

17 

0.978 

291.83 

0.10 

18 

0.995 

292.08 

0.02 
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Fig,  1  The  first  six  vertical  modes  of  the  18-level  NOGAPS  model, 
based  on  the  parameter  values  in  Table  1  (signs  are  arbitrary ). 


terized  physical  processes,  including  those  for  gravity- 
wave  drag,  vertical  fluxes  of  moisture,  heat  and 
momentum,  cumulus  convection,  large-scale  precipi¬ 
tation.  shallow  cumulus  mixing,  and  radiation.  Of 
particular  note  are  the  cumulus  convection  scheme, 
which  is  a  modified  version  of  the  Arakawa-Schubert 
(A-S)  scheme  (Arakawa  and  Schubert  1974;  Lord  et 
al.  1982).  and  the  vertical  flux  calculations,  which  are 
based  on  the  A7theory  formulation  of  Louis  ( 1979). 
Ground  temperature  is  predicted  over  land,  while  SST 
remains  fixed  at  its  initial  values.  The  latter  are  obtained 
from  daily  analyses  produced  by  Fleet  Numerical 
Oceanography  Center  ( FNOC ) . 

The  normal  modes  describe  the  spectrally  truncated 
mass  and  momentum  fields  and  are  derived  from  the 
adiabatic  model  equations  linearized  about  a  basic  state 
at  rest  with  a  horizontally  uniform  surface  pressure 
field  ps  and  a  horizontally  uniform,  vertically  varying, 
temperature  field  T( <r).  The  values  of p5  and  f(a)  used 
in  this  study  are  their  respective  global-mean  values 
averaged  over  days  21-50  of  a  control  simulation. 
There  are  18  vertical  modes  (corresponding  to  the 
number  of  vertical  levels)  numbered  /  =  1,  •  -  •,  18,  in 
in  which  /  =  1  denotes  the  gravest,  or  external  mode. 
The  calculated  values  of  f(<j)(or  7/)  used  in  this  study 
are  listed  in  Table  1 ,  while  the  vertical  structures  of  the 
six  gravest  modes  are  shown  in  Fig.  1 .  The  horizontal 
structures  of  the  modes  are  obtained  in  a  standard 
fashion  based  on  the  separability  of  the  linearized  sys¬ 
tem.  Further  details  are  given  in  Hogan  et  al.  ( 1991 ). 


3.  Numerical  experiments 

The  data  for  this  study  were  obtained  from  two  pairs 
of  50-day  simulations,  each  consisting  of  a  control  run 
and  an  anomaly  run.  The  two  runs  differed  only  in 
that  a  localized  positive  SST  anomaly  was  added  to 
the  initial  climatological  SST  values  of  the  anomaly 
run.  The  SST  anomaly  used  in  these  simulations  was 
located  in  the  western  Pacific  and  is  shown  schemati¬ 
cally  in  Fig.  2.  ft  is  equatorially  symmetric,  extending 
latitudinally  from  approximately  I4°N  to  14°S.  and 
longitudinally  from  105°  to  170°E.  The  magnitude  of 
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Fig.  2.  Schematic  representation  of  the  positive  temperature 
anomaly  added  to  the  climatological  SST  values  in  the  anomaly  run. 
The  outermost  contour  encloses  a  region  of  I  °C  anomalous  warming, 
while  the  innermost  contour  encloses  a  region  of  3°C  anomalous 
warming. 
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FiO.  3.  The  30-day  mean  12-h  accumulated  precipitation  at  1200  UTC  over  the  Pacific  Ocean  for  days  21-50  of  (a)  the  control  run  and 
(b)  the  response  field.  The  minimum  contour  in  (a)  is  I  mm  ( 12  h)*1,  with  contours  every  2  mm  ( 12  h)'1  thereafter.  In  (b)  the  contour 
interval  is  I  mm  ( 12  h)  with  negative  values  denoted  by  dashed  lines. 


the  anomaly  increases  from  +  1  °C  at  its  outer  boundary 
to  +3°C  at  its  center.  This  is  substantially  larger  than 
observed  SST  anomalies  in  this  region,  but  generally 
comparable  with  those  used  in  other  modeling  studies 
cited  in  section  1 .  The  inability  of  many  GCMs  to  show 
significant  responses  to  realistic  SST  anomalies  remains 
a  major  shortcoming  of  these  models,  but  is  not  of 


critical  importance  to  the  objectives  of  the  present 
study. 

The  full  SST  anomaly  was  introduced  at  the  begin¬ 
ning  of  the  simulations  and  kept  constant  in  time.  Al¬ 
though  there  should  be  no  significant  changes  in  the 
seasonally  averaged  model  response  as  to  whether  the 
anomaly  is  introduced  gradually  or  abruptly  some  time 
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prior  to  the  beginning  of  the  averaging  period  (e.g.. 
Mechoso  et  at.  1987),  the°'  may  be  some  differences 
during  the  early  stages  of  the  simulations.  In  particular, 
there  is  ev  idence  of  a  brief  period  of  model  shock  but 
no  abrupt  changes  in  the  rotational  component  of  the 
flow.  This  is  discussed  further  in  section  5. 

The  initial  conditions  for  the  control  run  in  the  first 
pair  of  simulations  were  based  on  the  NOGAPS  anal¬ 
ysis  from  1200  UTC  15  December  1985,  while  those 
for  the  second  pair  were  based  on  the  N'OGAPS  analysis 
from  1200  UTC  13  December  1987.  The  model  output 
was  archived  at  24-h  intervals,  so  that  the  output  data 
represent  the  instantaneous  values  of  the  fields  at  1 200 
UTC  on  each  day  of  the  simulation.  Composite  control 
and  anomaly  datasets  were  formed  by  averaging  the 
output  fora  given  number  of  days  from  both  the  1985 
and  1987  Northern  Hemisphere  wintertime  cases.  The 
results  that  will  be  discussed  in  the  following  sections 
are  based  on  these  composite  datasets. 

4.  Time-mean  response 

This  section  presents  some  of  the  general  character¬ 
istics  of  the  time-mean  model  response.  We  define  the 
response  as  the  difference  <5  between  the  control  C  and 
anomaly  .-I  simulations,  so  that.  8  =  .9  -  C.  The  time- 
mean  response  is  then  defined  as  the  composite  30- 
day  mean  value  of  8  for  days  2 1  -50  of  the  simulations. 
The  averaging  was  begun  at  day  21  in  order  to  ensure 
that  the  model  had  reached  an  approximate  equilib¬ 
rium  state.  Evidence  that  the  model  response  achieves 
equilibrium  by  day  20  is  presented  later. 

Figures  3a  and  3b  shows  the  30-day  mean  1 2-h  ac¬ 
cumulated  precipitation  from  0000  to  1 200  UTC  over 
the  Pacific  Ocean  for  the  control  run  and  the  response 
field,  respectively.  The  precipitation  pattern  for  the 
control  run  in  Fig.  3a  is  rather  typical  for  the  tropical 
Pacific  Ocean  during  the  Northern  Hemisphere  winter. 
There  is  a  well-defined  intertropical  convergence  zone 
(1TCZ)  centered  near  5°-10°N  that  spans  most  of  the 
Pacific  Ocean,  with  maximum  precipitation  rates  of 
7-10  mm  ( 12  h)  1 .  The  location  and  strength  of  this 
feature  are  reasonable  for  the  middle  of  winter  in  the 
Northern  Hemisphere.  There  are  also  locally  larger 
maxima  over  the  maritime  continent  and  near  the  ex¬ 
treme  west  coast  of  South  America.  The  precipitation 
rates  in  these  regions  are  on  the  order  of  15-20  mm 
( 12  h)  ' . 

The  response  field  shown  in  Fig.  3b  indicates  that 
precipitation  has  been  enhanced  substantially  in  the 
central  and  western  parts  of  the  region  of  increased 
SST,  especially  over  New  Guinea  and  the  extreme 
northeastern  part  of  Australia.  The  maximum  precip¬ 
itation  rates  in  this  area  range  from  5  to  10  mm  ( 12 
h ) ' 1  and  are  comparable  with  those  obtained  by  Kes- 
havamurty  ( 1982)  and  Blackmon  et  al.  ( l983)inGCM 
studies  of  SST  anomalies  in  the  central  Pacific.  These 


features,  combined  with  the  negative  precipitation 
anomalies  over  the  extreme  eastern  Pacific,  are  con¬ 
sistent  with  a  more  vigorous  Walker  circulation 

Figure  4  shows  the  area-averaged  30-day  mean  cu¬ 
mulus-heating  profiles  in  the  region  of  the  SS  I  anom¬ 
aly.  The  increased  SS  I  results  in  an  increase  in  cumulus 
heating  of  slightly  more  than  1  K  day  !  in  the  middle 
and  upper  levels.  The  dramatic  difference  in  the  char¬ 
acter  of  the  profiles  in  the  lowest  levels  is  due  primarily 
to  increased  precipitation  and  evaporative  cooling  in 
the  anomaly  run  compared  with  an  unrealistic  over¬ 
abundance  of  moisture  at  these  levels  in  the  con¬ 
trol  run. 

Figure  5  shows  the  30-day  mean  response  of 'he  ve¬ 
locity  potential  at  150  mb.  The  divergent  component 
of  the  wind  is  directed  perpendicuiar  to  the  contours 
toward  positive  values.  The  results  clearly  indicate  that 
there  is  increased  upper-level  horizontal  divergence  as¬ 
sociated  with  the  increased  convection  in  the  vicinity 
of  the  SST  anomaly.  This  pattern,  in  turn,  implies  in¬ 
creased  upward  motion  throughout  the  troposphere  in 
this  region,  with  convergence  and  downward  motion 
over  most  of  the  Atlantic  and  Europe.  This  is  consistent 
with  the  results  in  Figs.  3b  and  4. 

The  enhanced  Walker  circulation  is  also  evident  in 
Fig.  6.  which  shows  the  30-day  mean  vector  wind  re¬ 
sponse  at  1 50  and  925  mb.  In  Fig.  6a.  there  are  anom¬ 
alous  westerlies  at  150  mb  extending  across  most  of 


Fig.  4.  Area-averaged  30-day  mean  cumulus-heating  profiles 
over  the  region  of  increased  SST. 
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the  equatorial  Pacific,  indicating  flow  away  from  the 
SST  anomaly  lying  to  the  west.  These  anomalies  are 
on  the  order  of  10  m  s  '.  To  the  west  and  north  of 
New  Guinea,  there  are  easterlies  of  similar  magnitude. 
Note  the  anticyclonic  couplet  on  the  poleward  flanks 
of  the  SST  anomaly.  The  situation  is  reversed  at  925 
mb.  in  Fig.  6b,  where  the  flow  is  weakly  convergent 
over  New  Guinea  and  there  are  easterlies  on  the  order 
of  5  m  s  1  in  the  central  Pacific. 

Figure  7  shows  the  30-dav  mean  geopotential-height 
field  at  300  mb  for  the  composite  control  run.  Note 
that  the  results  exhibit  the  characteristic  features  of  a 
time-averaged  Northern  Hemisphere  wintertime  flow. 
In  particular,  the  zonal  asymmetries  in  the  Northern 
Hemisphere  height  field  indicate  that  the  major  mid¬ 
latitude  jets  occupy  their  usual  positions  along  the  east 
coasts  of  Asia  and  North  America.  As  expected,  the 
east  Asian  jet  is  stronger  and  has  greater  longitudinal 
extent  than  the  North  American  jet.  The  ridge  along 
the  west  coast  of  North  America  is  also  reasonably  po¬ 
sitioned.  In  contrast  with  the  Northern  Hemisphere, 
the  midlatitude  height  field  in  the  Southern  Hemi¬ 
sphere  exhibits  relatively  little  zonal  asymmetry,  owing 
to  the  lack  of  strong  differential  heating  from  land-sea 
temperature  contrasts,  especially  during  summer. 

Figures  8a  and  8b  show  the  30-day  mean  geopoten¬ 
tial-height  response  at  300  and  700  mb,  respectively. 
The  magnitude  of  the  response  is  comparable  with 
those  obtained  by  other  investigators,  although  the 
present  values  may  be  slightly  greater.  Presumably,  the 
response  might  be  weaker  if  we  average  over  many 


more  cases.  In  general,  the  largest  responses  occur  in 
the  Northern  Hemirphere.  This  is  consistent  with  the 
findings  of  Simmons  (1982)  and  Webster  (1982). 
among  others,  who  demonstrated  that  the  remote  re¬ 
sponse  is  usually  strongest  in  the  winter  hemisphere 
and  thus  probably  depends  on  properties  of  the  mean 
winds  in  the  subtropics. 

The  response  over  North  America  shows  a  familiar 
arcing  pattern  of  alternating  highs  and  lows.  This  pat¬ 
tern  (including  the  response  over  Europe)  is  very  sim¬ 
ilar  to  that  obtained  by  Shukla  and  Wallace  ( 1983). 
who  examined  the  response  of  a  climate  model  to  the 
composite  SST  anomaly  of  Rasmusson  and  Carpenter 
( 1982).  The  centers  of  the  features  in  Fig.  8  are  dis¬ 
placed  downstream  from  those  of  the  classical  PNA 
teleconnection  pattern  described  by  Wallace  and  Gut- 
zler  ( 1981 ).  In  this  case,  the  centers  are  located  close 
to  the  nodes  of  the  PNA  pattern.  Note  that,  in  general, 
the  largest  responses  occur  in  the  vicinity  of  North 
America. 

Figure  8  also  shows  that  the  midlatitude  response 
has  substantial  equivalent  barotropic  structure,  as  in¬ 
dicated  by  the  vertical  phasing  of  the  highs  and  lows 
and  their  increasing  amplitude  with  height.  This  struc¬ 
ture  difFers  from  that  of  the  tropical  response  in  Fig. 
6,  in  which  the  large-scale  flow  is  in  opposite  directions 
at  upper  and  lower  levels.  The  structure  of  the  response 
in  each  region  is  consistent  with  the  expected  predom¬ 
inance  of  either  shallow  or  deep  vertical  modes.  This 
is  examined  in  more  detail  in  the  remainder  of  this 
study. 


Fig.  6.  The  30-day  mean  vector  wind  response  over  the  Pacific  at  (a)  150  mb  and  (b)  925  mb.  The  maximum  vector  length 
corresponds  to  a  wind  speed  of  15  m  s  '  in  (a)  and  10  m  s  1  in  (b). 


5.  Normal-mode  analysis 

For  practical  purposes,  we  utilize  the  normal  modes 
derived  as  part  of  the  nonlinear  normal-mode  initial¬ 
ization  ( NNM  l }  scheme  in  the  NOGAPS  model.  Ac¬ 
cordingly,  we  should  note  that,  while  the  use  of  a  basic 
state  at  rest  greatly  simplifies  the  computation  of  the 


modes,  it  yields  vertical  structures  that  differ  somewhat 
from  those  obtained  with  a  more  realistic  ( strong,  ver¬ 
tically  sheared )  basic  current.  Therefore,  care  must  be 
taken  when  attempting  to  interpret  the  physical  pro¬ 
cesses  that  produce  a  particular  modal  response.  The 
use  of  more  realistic  basic  states  in  similar  normal¬ 
mode  decompositions  has  been  considered  by  other 
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investigators (Errico  1990,  personal  communication), 
but  unfortunately,  they  are  computationally  imprac¬ 
tical  and  are  beyond  the  scope  of  this  study. 

a.  Global  energy 

The  structures  of  the  modes  have  been  normalized 
such  that,  for  mode  J,  the  total  kinetic  plus  available 
potential  energy  Ej  is  given  by  the  simple  expression 

Ej  =  7  a  jaj ,  ( 1 ) 

in  which  ctj  is  the  amplitude  of  mode  J.  7  is  a  constant 
that  gives  Ej  dimensions  of  energy  per  mass,  and  an 
asterisk  denotes  a  complex  conjugate.  The  derivation 
of  ( I )  is  presented  in  the  Appendix  of  GIL  For  no w, 
it  is  of  interest  to  note  that  if  both  the  horizontal  and 
vertical  structure  functions  are  orthogonal,  then  the 
modes  contribute  independently  to  the  total  energy  E, 
and  the  quantity 

E=ZEj  (2) 

j 

is  equal  to  the  total  kinetic  plus  available  potential  en¬ 
ergy  of  all  the  modes.  The  vertical  modes  of  the  NO¬ 
GAPS  model  do  not,  in  general,  satisfy  an  appropriate 
orthogonality  condition,  and  thus,  E  +  E  in  this  study. 
It  can  be  shown,  however,  that  the  vertical  modes  do 
satisfy  this  condition  to  a  high  degree  for  all  but  the 
shallowest  equivalent  depths  so  that  E  is  a  close  ap¬ 
proximation  to  E  (e.g„  Errico  1989). 


b.  rime-varying  response  to  forcing 

To  obtain  the  response  energy  for  mode  J.  we  replace 
(  1 )  with 

Ej  -  y(Sirj)(So*  ).  (  3  ) 

where  oaj  =  iaj)A  -  t  «/)<  is  the  difference  in  the  am¬ 
plitude  of  mode  ./  between  the  anomaly  and  control 
simulations.  Note  that  E)  is  a  positive-definite  quantity 
not  to  be  confused  with  the  difference  of  the  energy 
between  the  control  and  anomaly  simulations.  This  is 
an  appropriate  diagnostic  quantity  in  the  present  con¬ 
text  because  it  represents  a  root-mean -square  weighted 
average  of  the  differences. 

As  an  overview  of  the  model  response.  Fig.  9  show  s 
the  total  energy  of  the  response  { i.e..  all  mode  types) 
as  a  function  of  vertical  mode  /  for  two  5-day  means 
corresponding  to  days  1-5  and  days  21-25  of  the  sim¬ 
ulations.  For  each  vertical  mode  /,  we  plot  the  quantity 

(4) 

Jl  X, 

where  the  summation  includes  only  modes  that  have 
equivalent  depth  ht  (denoted  by  the  set  St),  and  E*  is 
approximately  equal  to  the  total  energy  of  the  response 
defined  analogously  to  (2). 

For  future  reference,  note  that  the  set  «?/  does  not 
necessarily  include  all  modes  with  equivalent  depth  h, 
except  in  those  plots  depicting  the  total  ( i.e.,  rotational 
plus  gravitational)  energy.  The  quantity  E*  thus  rep- 


f 


resents  a  normalized  energy  that  describes  the  fractional 
contribution  of  a  subset  of  modes  to  the  total  energy 
of  the  response.  Also,  the  zonally  symmetric  modes 
have  been  excluded  from  these  figures  due  to  a  formal 
problem  in  defining  those  modes  corresponding  to  the 
mean  velocity  and  geopotential  fields  ( the  former  being 
zero  valued  in  the  linearized  model).  Although  meth¬ 


ods  exist  for  treating  these  modes  formally  (Kasahara 
1978;  Errico  1 987 ).  their  absence  has  little  appreciable 
impact  on  the  characteristics  of  the  response;  that  is. 
there  arc  only  small  differences  between  the  zonal-mean 
states  in  the  control  and  anomaly  simulations. 

Figure  9  shows  the  change  in  the  total  energy  dis¬ 
tribution  among  the  vertical  modes  between  days  1-5 
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and  days  21-25.  Duiing  days  1-5.  the  internal  modes 
dominate  the  response,  particularly  at  /  -  4.  as  well  as 
at  /  -  16  and  /  =  17.  By  days  21-25,  however,  the 
external  ( I  -  I )  modes  clearly  dominate  the  response, 
contributing  nearly  40r;  of  the  total  energy.  Only  those 
modes  corresponding  to  /  «:  5  contribute  significantly 
to  the  response  by  this  time,  with  the  strongest  medium- 
depth  response  shifting  from  /  -  4  to  /  =  5. 

The  responses  at  I  -  16  and  /  =  17  during  days  I- 
5  may  be  due  to  direct  thermal  forcing  by  extremely 
shallow  boundary -layer  processes  in  response  to  the 
introduction  of  the  SST  anomaly.  It  should  be  cau¬ 
tioned,  however,  that  the  available  potential  energy 
measure  that  has  been  derived  for  the  linearized  model 
is  not  the  same  as  for  a  nonlinear  mode!  (cf.  Lorenz 
1460b).  and  tends  to  exaggerate  the  significance  of 
modes  with  very  small  equivalent  depths  ( Errico  1987 ). 
In  any  case,  these  modes  are  of  no  physical  interest  in 
the  present  study. 

A  more  detailed  picture  of  the  evolution  of  the  re¬ 
sponse  for  vertical  modes  I  -5  is  shown  in  Fig.  10.  Here 
consecutive  5-day  mean  values  of  E]  are  plotted  for 
each  of  these  vertical  modes.  The  energy  contributions 
from  all  horizontal  mode  types  is  shown  in  Fig.  10a, 
while  the  contributions  from  the  rotational  and  grav¬ 
itational  modes  are  shown  separately  in  Figs.  10b  and 
10c.  respectively.  Note  the  different  ordinate  scaling 
for  the  gravitational  modes  in  Fig.  10c. 

Comparing  figs.  10a  and  1 0b.  we  observe  that  the 
strong  external  mode  response  is  comprised  almost  en¬ 
tirely  of  rotational  modes,  which  is  consistent  with 
idealized  modeling  results  suggesting  that  equivalent 
barotropic  ( and  primarily  nondivergent )  Rossby  waves 
play  an  important  role  in  the  response  ( Hoskins  and 
Karoly  1981;  Lim  and  Chang  1983;  Webster  1981). 
This  result  is  also  consistent  with  the  vertical  phasing 
of  the  geopotential-hcight  response  shown  in  Fig.  8. 


furthermore,  the  sinulariiv  of  the  growth  curves  atn-i 
the  first  week  in  f  igs.  10a  and  I  Oh  indicates  that  Oh- 
rotational  modes  dearly  dominate  th<-  response  at  ill 
vertical  scales  by  this  time. 

The  evolution  of  the  /  -  4  response  in  f  ig  10a  is 
quite  different — both  qualitatively  and  quantua 
lively — from  those  of  the  other  prominent  responses 
at  /  =  1  and  /  =  5.  The  /  =■  4  response  reaches  its 
maximum  value  during  days  1  -5  and  then  rapidly  ap¬ 
proaches  an  approximate  state  of  equilibrium  In  con¬ 
trast.  the  /  =  1  and  /  -  5  responses  reach  their  maximum 
values  two  to  three  weeks  into  the  simulation  \  com¬ 
parison  of  Figs.  1 0b  am)  10c  shows  that  the  /  4  re¬ 

sponse  in  Fig.  10a  is  primarily  a  gravitational-mode 
response  initially.  After  days  6-10.  however,  the  grav¬ 
itational-mode  energy  at/  4  decrease's !  Fig.  Ilk  i  while 
the  rotational-mode  energy  at  /  4  increases  toward 

its  equilibrium  value  (  Fig.  10b  i.  I  he  rapid  peak  in  the 
gravitational-mode  energy  is  undoubtedly  due.  in  pan. 
to  the  "shock"  effect  after  the  SSI  anomaly  is  intro¬ 
duced.  In  contrast,  the  responses  at  /  1  and  /  5  arc 

dominated  by  the  rotational  modes  throughout  most 
of  their  evolution. 

The  results  suggest  that  the  internal  responses  at  / 

4  and  /  5  may  be  of  different  dy  namical  origins 

Although  the  vertical  structures  of  both  modes  arc- 
consistent  with  the  "baioelinic"  type  of  response  as¬ 
sociated  with  tropical  convection,  in  which  the  large- 
scale  (low  reverses  sign  in  the  middle  troposphere,  the 
relatively  slower  growth  of/  5  and  the  predominance 
of  the  rotational  modes  throughout  its  evolution  suggest 
that  this  may  not  necessarily  be  the  case.  Furthermore, 
we  must  consider  that  the  nearly  barotropic  structure 
of  the  external  mode  in  the  troposphere  ( f  ig.  I  I  is.  by 
itself,  inadequate  for  resolving  the  equivalent  barotro¬ 
pic'  structure  of  the  midlatitude  response,  which  has 
increasing  amplitude  with  height.  Thus,  the  response 
at  I  -  5  ( maximum  amplitude  at  approximately  250 
mb )  may  instead  be  another  manifestation  of  the  mid- 
latitude  response.  This  possibility  is  examined  in  further 
detail  in  GIL  w  here  the  tropical  modal  response  is  sep¬ 
arated  from  the  midlatitude  one. 

The  results  in  f  ig.  10  demonstrate  that  a  significant 
global  response  occurs  rapidly,  with  all  vertical  modes 
shown  obtaining  ( at  least )  a  significant  fraction  of  their 
final  amplitudes  in  less  than  two  weeks.  An  important 
implication  of  these  results  is  that  the  time  scale  of  the 
global  response  to  tropical  forcing  is  clearly  relevant 
to  deterministic  NWP.  In  particular,  the  external  ro¬ 
tational  modes,  which  dominate  the  global  response, 
contribute  nearly  20'7  of  the  response  energy  after  only 
approximately  1  week,  and  30%-40*3  after  2-3  weeks. 
The  medium-depth,  internal  rotational  modes  also  play 
a  major  role,  especially  /  ~  5,  which  contributes  1  (Kv 
of  the  response  energy  after  approximately  one  week 
and  nearly  20'v  after  two  weeks. 
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6.  Horizontal  structure 

Further  insights  into  the  mechanisms  that  produce 
global  teleconnections  may  he  gained  by  examining 
the  dominant  horizontal  scales  of  the  model  response. 
Figure  1 1  shows  the  30-day  mean  300-mb  stream- 
function  response,  along  with  the  contributions  to  the 
response  from  those  modes  corresponding  to  zonal 
wavenumbers  1-3  and  5-9.  For  convenience,  these 
wave  groups  will  be  referred  to  as  the  long-  and  short¬ 
wave  components  of  the  response,  respectively.  Al¬ 
though  there  is  evidence  that  other  diagnostic  ap¬ 
proaches.  particularly  ones  that  account  for  two-di¬ 
mensional  wave  dispersion,  may  be  more  appropriate 
than  wavenumber  separation  for  examining  tropo¬ 
spheric  variability  (  Wallace  and  Hsu  1983 ).  the  short¬ 
comings  of  this  technique  are  somewhat  mitigated  by 
focusing  on  wavenumber  groups  rather  than  individual 
wavenumbers. 

F-igures  I  lb  and  1  Ic  reveal  a  dramatic  difference 
between  the  response  patterns  at  each  scale.  The  long 
waves  in  Fig.  I  lb  are  zonally  elongated  and,  in  some 
locations,  are  aligned  in  patterns  that  suggest  energy 


propagation  from  the  tropics  toward  ..igher  latitudes. 
The  most  prominent  of  these  patterns  originates  in  the 
far  eastern  equatorial  Pacific  near  Ecuador  and  appears 
as  a  wave  train  that  arcs  northeastward  along  the  east 
coast  of  the  United  States,  and  then  toward  southern 
Greenland.  The  amplitude  and  spatial  scale  of  the 
waves  increase  with  increasing  latitude,  consistent  with 
simple  theory.  There  are  hints  of  a  corresponding,  but 
noticeably  weaker,  branch  of  the  wave  train  in  the 
Southern  Hemisphere  that  arcs  southeastward  along 
the  west  coast  of  South  America.  This  difference  in 
amplitude  is  consistent  with  the  results  in  Fig.  8.  in 
which  a  generally  weaker  summer  hemisphere  response 
is  shown. 

The  large  longitudinal  separation  between  the  east¬ 
ern  Pacinc  wave  train  and  the  SST  anomaly  (Fig.  2) 
is  striking  and  raises  some  concern  as  to  whether  this 
may  be  only  a  stochastic  feature  of  a  relatively  small 
ensemble  average.  Webster  and  Chang  ( 1988),  how¬ 
ever.  have  demonstrated  that  longitudinal  trapping  of 
equatorial  wave  energy  can  produce  preferred,  geo¬ 
graphically  fixed,  regions  along  the  equator  for  pro¬ 
ducing  meridionally  propagating  responses.  A  corner- 
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FK(.  II  The  30-day  mean  3(H)- mb  stream  function  response  for  (a)  all  zonal  wavenumbers,  ( b )  wavenumbers  1-3.  and  (c>  wavenumbers 
5-9.  The  contour  interval  is  1. 5  x  l()6  m‘  s  1  in  fa)  and  1.25  x  10*  m?  s  1  in  (b)  and  fc),  with  negative  values  denoted  by  dashed  lines 


genceof  energy  occurs  in  these  locations  due  to  differ-  mary  location  for  the  accumulation  and  subsequent 
ential  changes  in  the  group  characteristics  of  tropical  meridional  propagation  of  equatorial  transient  energy 
modes  as  they  encounter  zonal  gradients  in  the  hack-  originating  from  as  far  away  as  the  western  Pacific.  The 
ground  flow.  These  authors  identified  the  eastern  Pa-  location  of  the  eastern  Pacific  wave  train  in  Fig.  I  lb  is 
ciftc  near  90°W — to  the  east  of  the  westerly  maximum  generally  consistent  with  this  scenario,  based  on  a  pre- 
in  upper-tropospheric  equatorial  winds — as  the  pri-  liminary  analysis  of  the  wind  field.  Moreover,  in  the 


present  case  there  is  also  a  substantial  positive  precip¬ 
itation  ( heating)  anomaly  in  the  eastern-centra)  equa¬ 
torial  Pacific  near  145°W  (Fig.  3b)  that  may  be  an 
effective  energy  source  for  this  response.  A  more  de¬ 
tailed  analysis  is  required  to  confirm  these  speculations. 

In  the  western  Pacific,  there  is  some  indication  that 
heating  in  the  vicinity  of  the  SST  anomaly  produces  a 
meridionally  propagating  response  in  the  Southern 
Hemisphere  that  emanates  from  just  west  of  Indonesia 
and  arcs  along  the  west  coast  of  Australia.  In  this  case, 
the  northern  branch  of  the  wave  train  is  conspicuously 
absent. 

In  contrast  to  the  long  waves,  the  short  waves  in  Fig. 

1  Ic  are  elongated  meridionally  rather  than  zonally  and 
are  aligned  in  a  zonal  pattern  between  latitudes  30° 
and  50°  in  both  hemispheres.  This  component  of  the 
response  bears  a  strong  resemblance  to  the  total  re¬ 
sponse  in  Fig.  1  la.  In  both  hemispheres,  there  are  zonal 
asymmetries  in  the  magnitude  of  the  response.  In  the 
Northern  Hemisphere  in  particular,  there  is  a  pro¬ 
nounced  increase  in  the  magnitudes  of  the  eddies  as 
we  move  across  the  Atlantic.  The  response  reaches  its 
maximum  amplitude  over  Europe  where  the  eddies 
have  more  than  twice  the  amplitude  of  those  in  the 
Pacific. 

Figure  12  shows  the  evolution  of  the  short-wave  re¬ 
sponse  at  days  7,  9.  and  10  of  the  simulation.  At  day 
7  (Fig.  12a).  the  first  organized  pattern  develops  with 
eddies  originating  off  the  California  coast  and  stretching 
across  the  United  States  and  into  the  eastern  Atlantic. 


The  rest  of  the  Northern  Hemisphere  is  relatively  un¬ 
disturbed  at  this  time.  In  the  Southern  Hemisphere, 
the  eddies  appear  to  originate  in  two  locations:  one  in 
the  eastern  Pacific  near  1 10°\V,  and  the  other  in  the 
extreme  southwestern  Indian  Ocean.  By  day  9  (Fig. 
12b),  the  response  in  the  Northern  Hemisphere  has 
propagated  eastward  and  equatorward  into  northern 
Africa.  The  eddies  in  the  Northern  Hemisphere  con¬ 
tinue  to  be  strongest  as  they  move  into  the  Atlantic. 
In  contrast,  the  response  in  the  Southern  Hemisphere 
attains  its  full  longitudinal  extent  more  rapidly  due,  at 
least  in  part,  to  the  existence  of  the  two  apparent  source 
points.  By  day  10  (Fig.  12c).  the  response  in  the 
Northern  Hemisphere  has  spread  across  the  Arabian 
peninsula  and  up  into  Asia,  leaving  only  the  eastern 
Pacific  relatively  undisturbed  The  response  develops 
on  a  time  scale  that  is  consistent  with  the  ones  observed 
in  Fig.  10.  which  shows  that  the  dominant  modes  grow 
rapidly  during  the  first  two  weeks  of  the  simulation. 

Although  the  long-  and  short-wave  partitionings  used 
here  are  somewhat  arbitrary  ,  the  results  exhibit  features 
that  are  consistent  with  some  current  theories  of  at¬ 
mospheric  teleconnections.  In  particular,  they  suggest 
that  the  background  flow  may  play  an  important  role 
in  augmenting  the  extratropicai  response.  This  is  es¬ 
pecially  compelling  in  the  Northern  Hemisphere,  where 
the  eddies  develop  near  the  exit  region  of  the  strong 
Pacific  jet  (cf.  Fig.  7).  Recently.  Murphree  et  al.  ( 1991 ) 
have  shown  that  barotropically  unstable  regions  on  the 
equatorward  flank  of  the  subtropical  or  midlatitude  jet 


Fig.  12.  The  3(J0-mb  Mreamfunction  response  for  zonal  wavenumbers  5-9  corresponding  to  (a)  day  7. 
( b )  day  9.  and  ( c )  day  1 0  of  the  simulations.  The  format  is  the  same  as  in  Fig.  1 1  c. 


can  be  the  dominant  energy  source  for  producing  an 
extratropical  response  to  a  tropical  heating  anomaly. 
Using  the  wave-energy  flux  vector  developed  by  Plumb 
(1985),  these  authors  showed  that  poleward-propa¬ 
gating  energy  from  a  tropical  heat  source  may  act  pri¬ 
marily  as  a  trigger  mechanism  for  this  instability. 


which,  in  turn,  produces  the  type  of  wave  response 
associated  with  midlatitude  atmospheric  teleconnec¬ 
tions. 

In  the  present  case,  it  is  also  interesting  to  note  that 
the  eddies  in  Fig.  I  lc  obtain  local  maxima  just  down¬ 
stream  from  where  the  meridionally  propagating  re- 
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spouses  in  Fig.  1  lb  intersect  the  midlatitude  westerly 
flow.  Conversely,  note  that  the  absence  of  a  local  max¬ 
imum  in  the  North  Pacific  short-wave  response  is  con¬ 
sistent  with  the  absence  of  a  northern  branch  of  the 
meridionally  propagating  response  emanating  from 
Indonesia. 

7.  Conclusions 

Rapid  global-scale  interactions  between  the  tropics 
and  extratropics  were  investigated  by  analyzing  the  re¬ 
sponse  of  the  U.S.  Navy’s  operational  forecast  model 
to  localized  SST  anomalies  in  the  tropical  Pacific.  The 
model  normal  modes  were  utilized,  in  addition  to  con¬ 
ventional  difference-field  diagnostics,  to  analyze  the 
structure  and  evolution  of  the  dominant  horizontal  and 
vertical  modes  of  response.  By  monitoring  the  energy 
growth  in  these  modes,  it  was  shown  that  the  qualitative 
and  quantitative  character  of  the  long-term  response 
is.  to  a  large  extent,  well  established  after  only  one  to 
two  weeks  of  simulation.  The  results  indicate  that  the 
global  atmospheric  response  to  tropical  forcing  can  be 
rapid  and  may  involve  a  series  of  dynamic  mechanisms, 
especially  in  the  extratropics. 

The  normal-mode  analysis  revealed  that  the  external 
rotational  modes  dominated  the  global  response  after 
the  first  few  days,  with  substantial  secondary  responses 
at  vertical  modes  4  and  5.  Once  equilibrium  was  es¬ 
tablished,  the  external  modes  accounted  for  over  30% 
of  the  response  energy,  while  vertical  modes  4  and  5 


accounted  for  approximately  10^  and  18%,  respec¬ 
tively.  The  strong  contribution  from  the  external  modes 
is  consistent  with  the  predominantly  equivalent  baro- 
tropic  extratropical  response,  as  demonstrated,  for  ex¬ 
ample,  by  the  vertical  structure  of  the  30-day  mean 
geopotential-height  response.  It  should  be  noted,  how¬ 
ever,  that  the  external  modes  used  in  these  calculations 
have  predominantly  barotropic  ( rather  than  equivalent 
barotropic)  vertical  structure  because  of  the  simple  ba¬ 
sic  state  used  in  their  derivation.  Thus,  they  are  inad¬ 
equate  for  describing  phenomena  that  are  strongly 
equivalent  barotropic.  This  undoubtedly  accounts  for 
some  portion  of  the  response  by  the  medium-depth 
internal  rotational  modes,  which  have  their  maximum 
amplitudes  high  in  the  troposphere.  This  conclusion  is 
consistent  with  the  growth  rates  for  these  modes  being 
similar  to  those  of  the  external  modes,  especially  for 
vertical  mode  5.  In  GII.  a  technique  is  presented  for 
separating  the  tropical  and  extratropical  modal  re¬ 
sponses  in  order  to  address  some  of  these  issues. 

A  decomposition  of  the  extratropical  response  in 
terms  of  zonal  wavenumber  groups  revealed  two  dis¬ 
tinctly  different  types  of  wave  patterns.  The  first  type, 
which  corresponded  to  wavenumbers  1-3,  was  char¬ 
acterized  by  zonally  elongated  waves  aligned  in  patterns 
suggesting  distinct  regions  of  meridional  energy  prop¬ 
agation  from  the  tropics  toward  higher  latitudes.  In¬ 
terestingly,  the  most  prominent  of  these  patterns  was 
separated  longitudinally  from  the  strongest  tropical 
heating  response.  Certain  aspects  of  these  results  are 
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consistent  with  the  two-stage  teleconnection  scenario 
proposed  by  Webster  and  Chang  (  1988).  in  which 
transient  tropical  energy  can  propagate  along  the 
equator  and  become  longitudinally  trapped,  thereby 
creating  favorable  emanation  points  for  a  meridionally 
propagating  response. 

The  second  type  of  wave  pattern,  which  correspon¬ 
ded  to  wavenumbers  5-9  and  bore  a  strong  resem¬ 
blance  to  the  total  response,  gave  further  insights  into 
mechanisms  that  may  be  important  for  producing  the 
extratropical  response.  The  growth  rates  of  these  waves 
were  consistent  with  those  of  the  dominant  normal 
modes,  but  showed  a  wave  pattern  that  developed  in 
situ  and  propagated  more  or  less  zonallv  in  the  mid¬ 
latitudes.  The  evolution  of  the  pattern  suggests  that  the 
background  flow  may  play  an  important  role  in  pro¬ 
ducing  an  extratropical  response  to  tropical  heating 
anomalies,  as  suggested  by  other  investigators,  includ¬ 
ing  Simmons  et  al.  ( 1983)  and.  more  recently,  Mur- 
phree  et  al.  ( 1991 ). 

The  result!)  reported  here  raise  many  questions  that 
cannot  be  answered  satisfactorily  without  further  in¬ 
vestigation.  In  particular,  it  would  be  useful  to  examine 
a  larger  number  of  cases  in  order  to  establish  the  ro¬ 
bustness  of  features  revealed  by  the  wavenumber  de¬ 
composition.  Similarly,  no  attempt  was  made  to  ex¬ 
amine  the  role  of  the  mean  winds  in  the  tropics  and 
subtropics  in  controlling  the  locations  of  the  different 
wave  responses,  or  their  dependence  on  seasonal  vari¬ 
ations.  A  more  rigorous  investigation  of  how  the  ex¬ 
tratropical  response  depends  on  interactions  with  the 
background  flow  is  planned  for  future  work.  Also,  the 
simple  basic  state  used  to  derive  the  normal  modes  can 
be  a  complicating  factor  when  trying  to  interpret  how 
certain  physical  processes  project  onto  the  modes.  The 
use  of  a  more  realistic  basic  state  might  simplify  this 
problem,  but  would  make  computing  the  modes  pro¬ 
hibitive. 

Notwithstanding  the  above,  these  results  clearly 
demonstrate  the  complexity  of  global-scale  interactions 
resulting  from  tropical  forcing  and  the  rapid  impact 
of  these  interactions  on  the  general  circulation  over  a 
broad  range  of  spatial  and  temporal  scales.  With  regard 
to  NWP,  the  results  imply  that  the  prospect  for  ex¬ 
tending  our  deterministic  forecast  capability  beyond 
the  medium  range  depends  largely  on  our  ability  to 
eliminate,  or  at  least  slow  the  growth  of,  systematic 
errors  in  the  tropics.  This  conclusion  is  similar  to  ones 
made  by  Tiedtke  ( 1984)  and  Sardeshmukh  and  Hos¬ 
kins  ( 1988),  who  determined  that  uncertainties  asso¬ 
ciated  with  present  convective  parameterization 
schemes  may  limit  our  prospects  for  making  skillful 
long-range  numerical  forecasts.  Finally,  these  results 
demonstrate  the  importance  of  developing  new  tools 
to  diagnose  the  behavior  of  complex  environmental 
models.  This  is  the  topic  of  Part  11  of  this  study  (GII), 


in  which  we  extend  the  normal-mode  analysis  pre¬ 
sented  here  in  order  to  investigate  the  structure  and 
dynamics  of  the  tropical  and  extratropical  responses 
in  more  detail. 
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